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SUMMARY 


The Prandtl-Glauert method lor subsonic potential- 
flow of a compressible fluid has generally been believed 
to lead to an increase in the press ures ov er a slender 

body of revolution by a factor 1 jsj 1 - M]_^ (where Mi 
is Mach number in undisturbed flow) as compared with the 
pressures in incompressible flow. Recent German work on 
thi s probl em has indicated, however, that the factor 

l/^l - Mp2 is not applicable in this case. In the 

present discussion a more careful application of the 
Prandtl-Glauert method to three-dimensional flow gives 
the following results: 

The Prandtl-Glauert method does not lead to a 
universal velocity or pressure correction formula that i a 

independent of the shape of the body. The factor i/VT - Mi 
is applicable only to the case of two-dimensional flow. 

The increase with Mach number of the pressures over 
a slender body of revolution Is much less rapid than for 
a two-dimensional airfoil. An approximate formula from 
which the increase can be estimated is derived theoreti- 
cally. 

The increase with Mach number of the maximum axial 
interference velocity on a slender body of revolution in 
a closed wind tunnel Is given approximately by the factor 

1 /Cl - Mi 2 )^ 2 , rather than by the factor 1 /(l - Mi 2 ) 2 
previously obtained by Goldstein and Young and by Tsien 
and Lees. 
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INTRODUCTION 


Because of its simplicity and applicability to both 
two- and three-dimensional flow, the Prandtl-Glauert 
method has been utilized to obtain approximate solutions 
for a wide variety of problems in; subsonic compressible 
flow. The effects of wind-tunnel interference on thin 
airfoils and on slender bodies of : revolution in subsonic 
flow have been thoroughly discussed in references 1 and 2. 
In common with many earlier writers on the subject, how- 
ever, including Prandtl and von Karman (see references 3 
to 6), the authors of references 1 and 2 state without 
clear proof that the Prandtl -Glauer t method leads to an .■ 

Increase by a factor of 1 /\jl - M in the pressures 
acting on the surface of a slender body of- revolution . in 
an unbounded, uniform parallel flow. (The symbol Mi 
denotes Mach number In undisturbed flow.) 

The present paper was prompted by recent German work 
on the problem of the subsonic flow over a body of revolu- 
tion, which indicates that application of the factor 

1 /\fl j is entirely incorrect in this case (refer- 

ence a\ 7 and 8). Some German writers (references 9 and 10) 
have even stated that there is no effect of compressibility 
on the pressures over a slender body of revolution in an 
unbounded stream - at least in firstr^appfbximation. In 
the present paper, the general problem -of the nearly 
parallel . subsoni c comnresslble flo.w over a closed body 
is treated as a straight-forward boundary-value problem. 

It is pointed out that previous writers (references b 
to 6) were lad to erroneous conclusions for the case of 
three-dimensional flow either because they disregarded 
the boundary conditions or because they did not examine 
the boundary conditions with sufficient care. In order 
to illustrate the general method of the present paper, 
the problem of the subsonic flow over a slender body of' 
r.evo-lution is discussed and the particular case of the 
ellipsoid of revolution is studiedin detail. 

The question of the apoli cation of the prandt 1-Glauert 
method to three -dimensional subsonic flow was .treated 
correctly for the first time by Gothert. In reference 7 > 
Gothert shows that the pressure coefficients on the surface 
of a body in a nearly parallel subsonic flow are obtained 
by calculating the incompressible flow over a body, the 
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lateral coordinates of which are contracted In the ratio 

U - Ml 2 : 1 as compared with the lateral coordinates of 
the original body, and then multiplying the pressure coef- 
ficients on the surface of this distorted body by the 

factor — -• The results obtained by Gothert's proce- 

1 - M i 2 

dure and by the method of the present paper will be equiva- 
lent in general. It is hoped that the reasons for the 
essential difference between the two- and three-dimensional 
flow will be brought out more clearly by the present 
treatment than by Gothert's procedure. 

The author is indebted to Dr. S. Katzoff for valuable 
discussions of this problem and criticism of the present 
paper . 


ANALYSIS 

Prandtl-Glauert Method for Subsonic Compressible Flow 


The prlncioal assumptions of the Prandtl-Glauert' 
method for subsonic compressible flow are as follows: 

(1) The Influence of the viscosity and conductivity, 
of the gas Is neglected. 

(2) The pressure is a continuous function of the 
density. 

( 3 ) Changes in pressure and density In the flow are 
small compared with the mean pressure and mean density: 
that is, the deviations from a uniform parallel flow are 
small . 

The Prandtl-Glauert method thus deals with subsonic 
flows that are shock free, or isentropic, and therefore 
Irrotational by Helmholtz' laws of vortex motion for a 
gas. In a Cartesian coordinate (x, y, z) system, the 
velocity potential $ is expressed by the relation ■ 


¥ 


§(x, y, z) = U]_x + cp(x, y, z) 


( 1 ) 


* 


5 



NACA TN No. 1127 


where cp is the disturbance potential and Ui is the 
velocity of the undisturbed flow.' The corresponding 
velocity components are 


v x = I T i + u, v 7 = v, v z = w 


( 2 ) 


where 


og> da 6 cp 

U = — , V = — •*-> W = — - 

c-x dy dz 


( 3 ) 


and 


JH. JL_ 

71 > 77 > 71 <C\ 1 « ’ J 

Ui u x Uj 


With the assumption that the squar.es and products of the 
disturbance velocities u, v, w are negligible compared 
with the velocities themselves, the exact equation of 
continuity is approximated by the partial differential 
equation for the disturbance potential 



where 


Mi Mach number in undisturbed 1 flow — 

\ a l 

and 

a^ speed of sound in undisturbed flow 

When the Mach number is vanishingly small, equation ( 4 .) 
reduces to Laplace's equation for incompressible potential 
flow . 
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A solution of equation ( I4. } must satisfy the following 
boundary conditions: 

( 1 ) The disturbance velocities u, v, and w must 
vanish in the undisturbed flow. 

( 2 ) At the surface of a solid boundary, the component 

of the gas velocity normal to the surface must vanish. 
Since + u ~ Uq, this condition is 


Uq cos a + v cos p + w cos y = 0 (5) 


where a, p, and y . are the angles that the normal to 
the body surface makes with the x-, y-, and z-axea, 
respective ly . 

In order to solve this boundary- value problem, 
Prandtl and Glauert introduce new Cartesian coordinates 
c, rp £ which are related to the coordinates x, y, and z 
by the affine transformation 


£ = x 


■n - y 



( 6 ) 


£ = z \jl - M]_ 2 



In the new 9, rj, £ system," the disturbance potential; 
satisfies Laplace's equation 


dfa + + 6^2. _ ■ 

dg 2 6-n 2 . 6 t Z 



If the velocity potential is regarded as unchanged 
in the transformation, the value of cp at the point 
(x, y, z) in physical space is identical with the value 
of cp at the point (9, rj, £) in the new space given 
by equations (0) . The relations between the disturbance 
velocities in the £, i], £ system (u 1 , v 1 , w') and 
the physical disturbance velocities 'are as follows; 
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u(x. 

7 } 

z) 

= §£ = 
6x 

6cp 

= 

& 

dx 

v(x. 

7 > 

z) 

_ . 

_ 

dp 


t>7 

dp 

dy 

w(x. 

7, 

z) 

~ 2$ - 

6cp 


6z 


dz 


u»(£, ri, t ) 


Vi - Ml 2 v(|, n, £) l ( 3 ) 


\/l - Mi 2 w'(g, tj , i) 


The lateral coordinates in the £, p, £ system are con- 
tracted, but the lateral velocity components in the 
g, p, t space are magnified as compared with the lateral 
velocity components at the corresponding point— in physical 
space . 


The boundary conditions that must be satisfied by 
the solution of Laplace's equation (equation (7)) can now 
be formulated as follows: 


(1) In the undisturbed stream, u», v', and. w 1 must 
vanish. 


(2) At a point on the surface of the solid boundary 
(x 3 , y s , z 3 ), the lateral velocity components v s and w a 
must satisfy equation (5)* These velocities are given by 



where 


= * 


s > 


Li 


= Vi - y s . £ b = \T- “I 2 z - 
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The potential cp(g, r), t) is set up as a general solution 
of Laplace’s equation that satisfies condition (1) and is 
appropriate to the problem. The arbitrary constants in 
the general solution cp(g, r\, £) are determined by 
substituting equations (9) in equation (5). Thus, by use 
of the Prandtl-Glauert method the problem of the nearly 
parallel subsonic compressible flow over a closed body is 
reduced to a problem in incompressible potential flow. 

The procedure outlined in the present paper was 
strictly carried out for the first time by Bilh&rz and 
Holder in 19l|0 for the particular case of a body of 
revolution of 10 -percent thickness ratio, with its maximum 
thickness at I 4.0 percent chord (reference 11). Bilharz 
and Holder, however, did not discuss the general implica- 
tions of their results. Previous writers on the subject 
attempted to draw conclusions directly from the differen- 
tial equation ( Ip ) and overlooked the necessity of satis- 
fying the boundary conditions . 


Prandtl-Glauert Method for Two-Dimensional Subsonic Plow 


For two-dimensional subsonic flow, the Prandtl- 
Glauert method leads directly to the result that the 
pressures and the axial disturbance velocities on a thin 
cylindrical body at a given small angle of attack are' ~ 

increased by a factor of 1 /fyl - Mi 2 over the corre- 
sponding values for incompressible flow. The analysis 
in this case is simple largely because the contraction 
of the coordinates of the body surface, which is required 
by equations (9)* never enters the problem. For two- 
dimensional flow, the lateral velocity component on the 
surface of a thin cylindrical body v s (x s , y s ) differs 
from the velocity on the chord v(x s , 0 ) by a quantity 
of higher order. This fact can be seen from physical 
considerations or from the development of the velocity 
v(x a , y) in a Taylor’.s series’ 


v(x 


s > 


y) = y(x 


s >■ 


0 } + 



From the approximate continuity equation 
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and the fact that 

( ( " >u } 

V dx / 0 

is of the order of y a , 

it 


\ /Xg,0 



can be seen that 

v 3 (x s , y s ) 

differs from v(x 3 , 0) 

by 


a quantity of the order of v s 2. 

The boundary condition that the resultant velocity 
on the surface must be tangential to the surface can then 
be satisfied by the velocity on the chord v(x s , 0). This 
condition is simply 


v(x s , 0) = Uq tan (3 (x 3 ) 


( 10 ) 


where 



In equation (10) tan (3 (x s ) is the slope of the surface 
at the point (x 3 , y s ) and c - is the chord of the body. 
The boundary condition f or_ a solution of Laplace's equa- 
tion in the p plane is 

Uq tan (3 (£ s ) 

v'(§ 8 , 0) = 


where 

o < 

2 . = = . 2 _ 


\/l - Mi 2 
~ < c 


The appropriate solution of Laplace's equation is 
given by a continuous distribution of 11 singularities" 
(vortices and/or sources and 3inks) along the chord. The 
strength of these singularities pep unit length is directly 

proportional to the factor Uq/^/l - Mq2 in the expression 
for the required normal velocity v ' ( § s , 0) given by 

equation (11). In the limiting case of - an incompressible 
fluid. Mi approaches 0 and the’ strength of the singu- 
larities is oroportional to Uq . In the present case. In 
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which M]. 7^ 0, the vortex or the source-sink strength 
along the chord must therefore be greater than the 

strength for M x = 0 by a factor of l/^l - Mi 2 . Since . 
the axial disturbance velocity u*(g s , 0) is also pro- 
portional to the strength of the singularities, this 

velocity is also increased by the factor 
equation (8), therefore, 

ju(x 3 , 0)] Mi [u.(g a , 0)]^ , 

[u(x a , 0)], [u'(Ss, 0)] 0 ]j ± _ M] 2 

regardless of the shape of the body or the angle of 
attack ^ 

Integration of Bernoulli's equation gives 


l/\jl - Mi 2 . By 


where 


= 1 ' 
Cp G ]/l - M X 2 


p " p i 
1 \ 2 ' 
2 P 1 U 1 


and 

p local static pressure 

p^ static pressure in undisturbed flow 

p-|_ density in undisturbed flow 


( 13 ) 


9 
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Prandtl-Glauert Method for Three-Dimensional 

Subsonic Plow 


For three-dimensional subsonic flow, the application 
of the Prandtl-Glauert method is not quite so simple as 
for two-dimensional flow. The lateral velocity components 
on the surface of the body v a and w 3 do not differ 
from the lateral velocities on the axis of the body by 
quantities of higher order. In the case of a body of 
revolution, for example, the radial velocity component v r 

fix) 

near the axis is very nearly • equal to -r * — S where 

q.Trr 

f(x) represents the strength of the source-sink distri- 
bution along the axis and r is the distance from the 
axis. In fact, it can be shown from the equation of conti 
nuity (equation ( ip } ) and the condition of Irrotationality 
that the radial velocity component v r near the axis 
is of the form ^ 

V = £iiLl + g(x)r log r + ... 
r 

Thus, the error made by neglecting the additional terms 
In the expression for v r is In the ratio r2 log r 
f(x) 

to — and this error is beyond the limits of a first 

r 

approximation . 

The radial velocity component thus increases Indefi- 
nitely as the axis is approached -in contrast to the 
situation in the case of two-dimensional flow. Since the 
lateral velocity components on the surface of the body of 
revolution v 3 and w s obviously cannot be replaced by 
the lateral velocities on the axis, the contraction of 
the coordinates of the body in the p, l system enters 
the problem; consequently the shape of the body becomes 
important. This fact- has apparently been overlooked in 
the discussions of the problem given in references 1 to 6 . 

The boundary condition at a point on the surface of 
the body (x s , y s , z s ) is 
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v r s (x 3 , y s . 



- U]_ tan (3 


(x s ) 


(111-) 


where 


£ < 
2 ~ 


< £ 

- 2 


q on ( 1 J 4 .) tan (3 (x s ) la the slope of the surface 
a raeridan plane and c is the length of the bodv 
The boundary condition to be satisfied at the point 
vss* r)g, c, s J by a solution of Laplace's equation is 


v ’r g (£ s , n s , £ a ) = _. P (S »I 

]/l - Mi 2 


(15) 


where 


£ < x, < o 

2 - Ss ~ 2 


w e Q a J Dr ?? rlate aolutlon of Laplace's equation is given 
Y ontinuous distribution of sources and sinks along 

the axis between £ = -f and £ = of strength f(£) 

per unit length. The strength f(£) is directly related 
to the required lateral velocity v> ' by J reXated 


v, r s (£ s , Tis, t 3 ) = -L 


£(Zs) 




■= S 


f(§s) 


(16) 


Vi - Mx 2 fJTTJ 
^? ati ° n ( l 5) the s , l ength K£) Of the source. 

^nk QiscriDu tion o n cne axia is propor tional t o — 

1 ; l0r the glVen b 0 Q y/ inde pendently of the if a cl 

number or the undisturbed flow*! 
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(Since the equivalent source-sink distribution along 
the axis of a slender body of revolution is independent 
of Mach number in first approximation, an analysis similar 
to that given in reference 2 shows that the increase with 
Mach number of the maximum axial interference velocity on 
a slender body of -revolution in a closed wind tunnel is 

given approximately by the factor if (l - Mq^5/2. ra ther 
than by the factor 1 j(l - previously obtained in 

references 1 and 2. This result is in agreement with the 
conclusions reached by von Baranpff in reference 12.) 

The fact that f(£.) is independent of Mach number 
does not mean that the axial disturbance velocity 
u(x s , y s , z s ) = u'(£ 3 , p s , £ s ) is also independent of 

Mach number. The point (£ s , r) 3 , £ s ) lies closer to the 
axis for a Mach number Mq / 0 than for incompressible 
flow. For the given source- sink distribution, the value 
of u l (£s, rj a , L s ) near the maximum thickness of the 

body will therefore be larger than the corresponding value 
for incompressible flow. An estimate of this increase 
can be obtained by an examination of the incompressible 
field of flow around a thin ellipsoid of revolution 
(prolate spheroid). (See reference 13 . ) The disturbance 
velocity potential is . 

q>(X., p) = Ap^| log - l) ( 17 ) 


where \ and p are elliptic coordinates and A is a 
constant. (The symbol \ is used Instead of Lamb 1 3 £ 

in order— to avoid confusion.) in the vicinity of the 
maximum thickness of the body, : ■ 


\ — > 



+■ 


+ i 2 - 

(c/2) 2 


p — ^0 


> 


( 18 ) 


In equation (18) c is the length of ellipsoid of revolu- 
tion. The maximum axial disturbance velocity r| , £) 

is approximated by ; : . 
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uf max(£=* Tjj t) 


V 1 - p 2 / Ocp\ 

| - m- 2 W/^—^o 


~ 1 - log 2 + log 


V-n 2 + j c 

c/2 


( 19 ) 


Therefore 


u m.ax( x sj Ys 

where t/c is 
follows that 


, z s ) si- log 2 + log (Vi - Mi 2 
the thickness ratio of the body. It 


( 20 ) 


Wm! C PMi log v/l - Ml 2 

( u max) 0 C p ^ 0 . 31 + log | 


( 21 ) 


Schmieden and Kawalki in reference 8 obtained the same 
result by a careful application of the procedure outlined 
by Gothert (reference 7)* The original expression given 
by Gothert is Incorrect. 

Prom equation (21) the effect of compressibility can 
be seen to vanish as the thickness ratio approaches zero - \ 
In contrast to the case of two-dimensional flow - where V 
the first-o rder effect of compressibility is given by ' 

l/yl - Mq 2 , independently of the shape of the body or the ) 
angle of attack of the body. 

For slender bodies of revolution that are symmetrical 
or nearly symmetrical longitudinally, equation (21) should 
give some indication of the first-order effect of compressi- 
bility. This equation is plotted In figure 1 for various 
values of the thickness ratio t/c. The effect of com- 
pressibilit y is ev idently considerably smaller than the ____ 

factor 1 /\jl - Mq 2 .would predict. Eor an ellipsoid of 
revolution of 10-percent thickness ratio the estimated 
increase in maximum negative pressure coefficient is only 
25 percent at a Mach number of 0.80. The Increase in 

13 
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maximum pressure coefficient obtained in reference 11 by 
an approximate numerical procedure . for a body of revolu- 
tion of the same thickness ratio with its maximum thick- 
ness at I4.O perc ent cho rd is percent. By comparison, 

the factor l/^l - predicts an increase of 67 percent 

in the axial disturbance velocity and pressure coef fi cient . 

In passing, it shoul d be no ted that the inapplica- 
bility of the factor 1 /\jl - for three-dimensional 

subsonic flow is indicated by the results of bhe calcula- 
tion of the subsonic flow over a sphere by the Janzen- 
Rayleigh method of iteration (reference 14) . The first- 
order effect of compressibility on' the pressures must be 
less for a slender body of revolution than for a thick 
body, such as the sphere. Since the maximum pressure 
coefficient on the sphere at the critical Mach number r 

Ml = 0.573 is les3 than the factor l/^l - Mi^ would 
predict (fig. 1), this factor cannot possibly give the 
correct result for a slander body of revolution. 

Because the axial disturbance- velocity near the 
maximum thickness of a slender body of revolution is 
small (“O.OJOXJi in the example of reference 11) and the 
increase of this velocity with Mach number is also small, 
the assumptions of the Prandtl-Glapert method are more 
nearly satisfied at considerably higher subsonic Mach 
numbers in three-dimensional f_low than in two-dimensional 
flow. For the same reasons, the so-called critical Mach 
number for a slender body of revolution is quite close 
to unity. This conclusion is generally applicable to any 
three-dimensional body, the lateral dimensions of which 
are small compared with the dimensions in the direction 
of motion, for example, the highly swept-back thin wing. 


CONCLUDING REMARKS 


An analysis of the Prandtl-Glauert method for., sub- 
sonic compressible flow around a slender body- of revolu- 
tion led to the following conclusions: 

1 . The method does not give a universal velocity or 
pressure correction formula that is independent of the-. 

shape of the body. The factor i/yi - Mi 2 .(where Mi 

lit. 
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is Mach number in undisturbed flow) is applicable only 
to problems of two-dimensional flow. 

2. The increase with Mach number of the pressure 
coefficients and axial velocities on the surface of a 
slender body of revolution is much less rapid than for a 
two-dimensional airfoil. An approximate estimate of this 
Increase for slender bodies of revolution that are sym- 
metrical or nearly symmetrical longitudinally can be 
obtained from the expression • - 


i + log 

(%iax ) 0 0.31 + log | 


where 1 %^ is maximum axial disturbance velocity and 
t/c is the thickness ratio of the body. (This expression 
is derived from a consideration of a thin ellipsoid of 
revolution.) 

3 . The so-called critical Mach number for a slender 
body of revolution is quite close to unity. This con- 
clusion is generally applicable to any three-dimensional 
body, the lateral dimensions of which are small compared 
with the dimension in the direction of motion, for example 
the highly swept-back wing. 

I 4 .. The equivalent source-sink distribution along the 
axis of a slender body of revolution is independent of 
Mach number in first approximation. Therefore, the ■ 
increase with Mach number of the maximum" axial inter- 
ference velocity on a slender body of revolution in a 
closed wind tunnel is given approximately by the factor 

l/(l - Mj2) 5/2 rather than by the factor if (l - 

previously obtained by Goldstein and Young and by Tsien 
and Lees . 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , February 13 , 19^6 
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Figure L- Increase of maximum negative pressure coefficient with Mach number 

for a thin ellipsoid of revolution. 








